Shell model study of single-particle and collective structure in neutron-rich Cr 

isotopes 
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The structure of neutron-rich Cr isotopes is systematically investigated by using the spherical 
shell model. The calculations reproduce well the known energy levels for the even-even 52-62 q^. 
and odd- mass 53-59qj. jj^gjgj^ g^j^jj predict a lowering of excitation energies around neutron number 
A'' — 40. The calculated B{E2; 2^ Of) systematics shows a pronounced collectivity around 
A'' — 40; a similar characteristic behavior has been suggested for Zn and Ge isotopes. Causes for the 
sudden drop of the 9/2^ energy in ^^Cr and the appearance of very low Oj states around A'^ = 40 are 
discussed. We also predict a new band with strong collectivity built on the Oj state in the A*" = 40 
isotope ^*Cr. 

PACS numbers: 21.10.Dr, 21.60.Cs, 21.60.Jz, 21. 10. Re 



I. INTRODUCTION 

The neutron-rich /p-shell nuclei far from the valley of 
stability are of particular interest in recent experimen- 
tal and theoretical studies To provide a satisfactory 
description of these nuclei, the challenge remains to un- 
derstand what mechanisms cause changes in nuclear shell 
structure as neutron number increases in nuclear systems. 
Theoretical calculations have questioned the persistence 
of the traditional magic numbers, which have been known 
to exist in stable nuclei. For example, ^^Ni is expected 
to be a double-magic nucleus as a consequence of the 
neutron subshell gap separating the /p-shell and gg /2 or- 
bital at the neutron number = 40 [3, 0, |1| . The size 
of this subshell gap is characterized by the excitation 
across the /p-shell and 59/2 orbital [1, @]- In addition 
to this, the presence of the positive-parity 39/2 orbital 
above the negative-parity fp shell strongly hinders Ip-lh 
excitations. This A^ = 40 shell gap, however, disappears 
when protons are added to or subtracted from ^*Ni Q . It 
was discussed in Refs. [1, Q that in the early mean- field 
calculations, a distinct shell gap that exists in the TV = 
40 nucleus ^*Ni disappears when quadrupole correlations 
are taken into account. Even for ^^Ni, it does not show 
a pronounced irregularity in the two-neutron separation 
energy as expected for a typical double-magic nucleus. 
It has been suggested [13] that a small B{E2,0^ 2f ) 
value is not a strong evidence for the double-magic char- 
acter. We may thus conclude that the double-magicity 
nature in ^^Ni is still controversial and remains an open 
question. 

On the other hand, collective structure is predicted 
to develop in neutron-rich nuclei. Zn and Ge nuclei are 
known to exhibit rather strong collectivity around A'^ = 
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40 [Tl|, [I^l. Both unusually low excitation of the first 
excited Oj state and strong enhancement of B{E2; 0^ 
2^) near A^ = 40 indicate a dramatic structure change at 
this nucleon number. Recently, we have shown that 
this characteristic behavior can possibly be understood 
in terms of rapid increase in the gg /2 proton and neutron 
occupation. However, in order to explain the structure 
change and enhancement in B{E2), we needed effective 
quadrupole matrix elements and large effective charges. 
Shell model calculations ;(| have also predicted an excited 
band in ^^Ni based on the Oj state. 

For 60,62,64Qj.^ shell-model calculations in the fpg-sheW 
space with the ^®Ca core (i.e. the '^'^Ca core and eight 
fj/2 frozen neutrons) have shown that these nuclei 
are strongly deformed with large quadrupole moments 
and large B{E2; 2^ — > 0^) values. The recent report 
on some even-even Cr isotopes has indicated a lowering 
of the 2+ energy beyond A'' = 34. The monopole inter- 
action discussed by Otsuka et al. [l^, [l7| should affect 
differently the protons occupying the vr/y /2 and vf^ /2 or- 
bitals, and should have an impact on the neutron single- 
particle spectra involved. However, it has been demon- 
strated that the ug^ /2 orbital should not be ignored in the 
description of A^ > 34 isotopes. In fact, the GXPFIA in- 
teraction [isj in the /p-shell space cannot describe the 
energy levels in ^^Cr and ^°Cr [l^, and the calculated 
excitation energies are found much higher than the ex- 
perimental ones. Sorlin et al. [1| suggested an onset of 
a sizable deformation in ^°Cr. Recent observation [l5[, 
however, has questioned this interpretation in terms of 
deformation. Rather, the low-lying energy systematics 
of neutron-rich Cr isotopes has been discussed in terms 
of a softening of nuclear shape with increasing neutron 
number. 

The low-lying 9/2 j*" states in odd- mass neutron-rich Cr 
isotopes are obviously outside the /p-shell model space. 
The recent observation .19j in some odd-mass neutron- 
rich Cr isotopes reveals that the 9/2^ state energy drops 
down considerably with increasing neutron number. The 
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FIG. 1; Comparison of experimental energy levels with the 



KB3G and the GXPFIA calculations for 



^Cr. 



sharp decrease of the 9/2^ state energy is considered as 
a clear indication for the monopole interaction, and the 
attractive force between the 1 gg /2 and 2p^ /2 neutron or- 
bitals pushes down the 139/2 orbital. The mean-field cal- 
culations suggested the presence of prolately-deformed 
rotational bands built on the 59/2 states in ^^'^^Cr, while 
the long-lived gg/2 isomeric state observed in ^^Cr was 
considered to be consistent with an oblate deformation 
[20l . [2l| . All these indicate that different shapes may co- 
exist and a shape change may occur in ^^Cr. The abrupt 
change in structure of odd-mass neutron-rich Cr isotopes 
can be better explained by the appearance of the gg/2 
orbital at low energy. The observation of the 9/2f state 
at low-excitation energy undoubtedly demonstrates the 
necessity of including the 179/2 orbital in the calculation. 

When these neutron-rich nuclei rotate faster, it is re- 
quired also to involve the 99/2 orbital to discuss the 
physics. The study of high-spin states in the even-even 
56, 58, 60 (^j. isotopes is important to gain information about 
the evolution of single-particle and collective excitations. 
A comparison of the high-spin data with current shell 
model calculations indicates that there is strong moti- 
vation for extending the interactions to include the neu- 
tron gg/2 orbital [l5|. Thus, a consistent approach to the 
structure study in this mass region would suggest that 
the gg/2 orbital should drive toward collectivity in the 
neutron-rich Cr isotopes for N > 34. 

In Fig. [U we compare the experimental energy levels 
of ^^Cr and ^°Cr with the present calculation, and with 
calculations using the GXPFIA and KB3G interactions. 



It has recently been discussed [is'] that the GXPFIA in- 
teraction [Tsl is inadequate to describe the neutron-rich 
Cr isotopes beyond N = 34. For example, for ^^Cr their 
calculated energy levels of the low-lying negative-parity 
states lie much higher than what arc seen in the data, and 
the discrepancy is further amplified in higher spin yrast 
states. While their results for the < 34 Cr isotopes 
reproduce the data, the agreement with the even-even 
^°Cr data is poor for high-spin states. To understand the 
single-particle and collective structure in neutron-rich Cr 
isotopes, shell-model calculations in a full fpg-sheW space 
are highly desirable. However, conventional shell-model 
calculations in the full fpg-sheW space are not possible 
at present because of too huge dimensions in the config- 
uration space. For neutron-rich Cr isotopes, we have to 
restrict the model space to the fpg model space of "^^Ca, 
namely, using the ^"Ca core with eight /7/2 frozen neu- 
trons. In the N = AO region, the d^/2 orbital may need 
to be taken into account to develop quadrupole collectiv- 
ity [Hi Hi HI, [2l|. However, recent observations in '^"Zn 
and *^Ge jlll| indicate the persistence of the A^ = 50 
shell gap, implying that the d^/2 orbital may not affect 
much the collectivity of the neutron-rich Cr isotopes with 
N < 50. 

In the present work, we perform large-scale spherical 
shell model calculations using the pairing plus multipole 
forces with the monopole interaction included [l^, 12^ . 
The fpg model space comprises the O/7/2, IP3/2, O/5/2, 
lpi/2 active proton orbitals and O/7/2, IP3/2, O/5/2, lpi/2, 
O59/2 neutron orbitals with eight /7/2 frozen neutrons. 
Structure of neutron-rich Cr isotopes will be investigated 
in detail. In particular, the interplay between single- 
particle and collective properties will be discussed. 

The paper is arranged as follows. In Section |TT1 we 
outline our model. In Section [1111 we present the results 
from numerical calculations and discuss them for a long 
chain of Cr isotopes. Finally, conclusions are drawn in 
Section HVl 



II. THE SHELL MODEL 

We start with the following form of Hamiltonian, 
which consists of pairing plus multipole terms with the 
monopole interaction included: 

H = Hsp + Hpg + Hp.^ + Hqq + Hoo + ^JiT" + ^^mc 
a J=0,2 Mk 

M M 

~ ^° XI XI {ah) A J MOO [ah) 

a<b JM 

+ XX ^-cH) X A^jMTKiab)AjMTKiab), (1) 
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FIG. 2: (color online) Effective single-particle energies Cq for 
even-even Cr isotopes. (Upper panel) calculation from the 
shell- model Hamiltonian with the monopole interactions, and 
(lower panel) without the monopole interactions. 



where b in the third and fourth terms is the length param- 
eter of harmonic oscillator. We take the J = and J = 2 
interactions in the pairing channel, and the quadrupole- 
quadrupole {QQ) and octupole-octupole (OO) terms in 
the particle- hole channel [Ij, HH . The monopole interac- 
tion is divided into two parts, namely the average T = 
monopole field H^^^ and the monopole correction term 
i?mc- The Hamiltonian ([1]) is isospin invariant. 

Our calculation has be en p erformed by using the AN- 
TOINE shell-model code [1^. We consider the ''^Ca core 
as mentioned above, and employ the single-particle ener- 
gies e/7/2 = 0.0, ep3/2 = 2.0, epi/2 = 4.0, e/5/2 = 5.5 (aU 
in MeV) in the calculation. The relative single-particle 
energies of /7/2, P3/21 P1/2, and /5/2 are taken from the 
excitation energies of the low-lying negative-parity states 
in '^^Ca, and the energy between e/7/2 and £^3/2 is deter- 
mined from the excitation energy of 3/2^ state in ^^Ca. 
The single-particle energy £39/2 = 4.5 MeV is chosen so 
as to reproduce the 9/2^ energy level in ^"^Cr. It is not 
surprising that the gg /2 orbit lies below the /s /2 orbit be- 
cause the 9/2+ level is lower than the 5/2~ level in ^^Ca 
[3l| . This lowering of the /2 would be attributed to the 
attractive T = 1 monopole interaction VT^^„ 

We adopt the following interaction strengths for the 



pairing plus multipole forces 

go = 17.89M, 52 = 152.24M^/3, 

X2 = 228.36/A5/^ X3 = 485.1M2 (in MeV). (2) 

The average monopole force -ff Jj"" is neglected because 
it does not affect the low-lying excited states of neutron- 
rich Cr isotopes. For the monopole correction terms, we 
use 



k 



T=0 



(/7/2,P3/2) = 0.5, 
fc^f (/7/2,Pl/2) = 0.4, 



uT=0 



(/7/2: 15/2) 
(^3/2,39/2) 



-0.6, 
-1.0(3) 



The repulsive T^fr/2-i^P3/2 and T^fr/2-'^Pi/2 monopole in- 
teractions can give rise to the correct order of energy 
levels in odd-mass Cr isotopes, and also lower the first 
excited Oj state for the even-even Cr nuclei near = 28. 
From comparison of the low- lying levels of "^^Ca with 
those of "Ni, one sees clearly that the 5/2j^ state in ^^Ni 
requires a monopole interaction between nfj^2 and 1^/5/2. 
This interaction also plays an important role in lowering 
the first excited Oj state around iV = 40. The sharp de- 
crease of the low- lying 9/2^ state observed along the 
isotopic chain from '^^Cr to ^^Cr is considered as a clear 
indication for the monopole interaction, and the attrac- 
tion between the 1^^3/2 and h'gQ/2 orbitals pushes down 
the 9/2^ state. Having considered the above monopole 
interactions, we choose the single-particle energy of the 
gg/2 orbital so as to reproduce the 9/2^ energy in ^^Cr. 
For calculations of spectroscopic Q-moments and B{E2) 
values, we take the standard effective charges e^r = 1.5e 
for protons and = 0.5e for neutrons. 

To get insight of the characteristic behavior of these 
neutron-rich Cr isotopes, we further consider mean- field 
approximations for our shell model Hamiltonian. The so- 
obtained effective single-particles arc the dressed parti- 
cles that carry information on the interactions. In Fig. [21 
we present the calculated neutron effective single-particle 
energies (SPE) [£]. The upper (lower) panel in Fig. [5] 
shows the effective SPE in the shell-model Hamiltonian 
(III) with (without) the monopole interactions ([3]). In the 
upper panel, we can see drastic variations of the effective 
SPE for each orbitals. 

First, starting from A = 28, the effective SPE for 
the (7g/2 orbital decreases quickly with increasing neu- 
tron number, and then becomes degenerate with the P3/2 
orbital in the range from A^ = 32 to A^ = 40. The sud- 
den drop of the first excited 9/2^ in the odd-mass nuclei 
53-59(^j, (^ggg later in Figs. [5]and[S]) is related to this be- 
havior of 59/2, while the pairing correlation increases the 
excitation energy of 9/2^ by about 1 MeV. This rapid 
decrease of the effective SPE is attributed to the strong 
attractive monopole interaction between the i'P3/2 and 
h'gg/2 orbitals in Eq. 

Second, we can see a change of the SPE level spacings 
among the orbitals around A^ = 40. This change is also 
due to the monopole interactions in Eq. ([3]). The repul- 
sive 7r/7/2-i'P3/2 and Trf7/2-^Pi/2 monopole interactions 
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FIG. 3: Experimental and calculated energy levels of positive-parity states for even-even isotopes Cr. 



give a correct sequence of the energy levels in odd-mass 
Cr isotopes. The attractive tt 77/2-1^/5/2 monopole inter- 
action plays an important role in lowering the first ex- 
cited state around N — 40. As already mentioned be- 
fore, the attraction between the i'Pz/2 a-nd 1^39/2 orbitals 
pushes down the 9/2^ state, and provides an explanation 
for the sharp decrease of the low-lying 9/2 j*" state from 
^^Cr to ^^Cr. These monopole interactions cooperatively 
give rise to the reduction of the effective SPE spacing. 
The deformation that develops around A'^ = 40 is due 
to the quasi-degeneracy of the effective SPE spacing and 
the collectivity of the intruder vgg/2 orbital. 

III. RESULTS AND DISCUSSIONS 

A. Energy levels 

Let us first present the shell model results for even- 
even Cr isotopes. Figure [3] shows a comparison between 
calculated results for positive-parity states and the ex- 
perimental data [3, [3l| for the even-even Cr isotopes 
52-62(^j,^ As one can see, the theoretical results describe 
satisfactorily the experimental energy levels. The feature 



that the first excited 2'{ energy is steadily decreasing as 
neutron number approaching A^ = 40 is correctly repro- 
duced. One can also see a local increase of 2^ energy at 
s^Cr. 

The calculated first excited Oj states for 52-56(^j. gj^ow 
a decreasing trend with increasing neutron number, al- 
though the calculated Oj energy for ^^Cr lies too high 
as compared to experiment. This discrepancy may be 
caused by the assumption of the A^ = 28 shell closure in 
the present calculation. The neutron excitations of the 
A^ = 28 core, which are neglected in the present study, 
would lower the Oj state in ^^Cr. The situation is quite 
similar to the anomalous behavior of the Oj state in ^''Fe. 
In fact, the full /p-shell model calculation ^2] has shown 
that the Oj states in ^^Cr and ^"^Fe are the 2p-2h band- 
head. The present calculation predicts the Oj energy 
levels in 58-62Qj. which have not yet been observed. It is 
interesting to see a particular prediction that the calcu- 
lated Oj level jumps up at ^"Cr with A^ = 36 while the 2+ 
level in the same nucleus decreases. The attractive i'p^/2- 
vgg/2 monopole interaction pushes the 179/2 orbital down, 
causing the single-particle spacing between the (79/2 and 
the Pi/2 orbital to increase. 

Figure 2] shows a comparison between the calculated 
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FIG. 4: Experimental and calculated energy levels of 
negative-parity states for ^^'^*Cr. 



and experimental energy levels for the negative-parity 
states in ^^Cr and ^*Cr [l5|. In ^^Cr, the odd-spin lev- 
els of 7^ , 9^ , III , 15j^, and 17^ are reproduced 
very well, while the even-spin levels are predicted. The 
37 states in ^^Cr and '^^Cr are predicted to lie around 
2.5 MeV. The 3j~ state is considered as an Ip-lh excita- 
tion from the negative-parity orbital to the positive 
parity orbital gg/2 because the 3j~ energies are compara- 
ble to the single-particle spacing 2.5 MeV between these 
orbitals. In addition, it is interesting to see a predicted 
formation of spin doublet states in both nuclei, as for 
instance the pairs lOj^-llj^, 12j^-13]^, lA^-W^, and IQi- 
17 1 . In ^*Cr, the even-spin levels of 6~ and 8~ are re- 
produced well, although the calculated odd-spin and 
7i levels are lower than the experimental ones. 

In Fig. [51 calculated results for the odd-mass Cr iso- 
topes 53-59 (^j, g^pg shown, and compared with experimen- 
tal data [1^, [3l[. Overall, the calculation reproduces the 
data well. In particular, the agreement for ss.s^q^. .^yj^}^ 
both negative- and positive-parity states is excellent. In 
^^Cr, most energy levels are described well; only the low- 
lying 5/2^ and 7/2^ states are calculated higher than 
the experimental ones. In these odd-mass Cr isotopes, 
the most striking feature is that the first excited 9/2f 
state drops down rapidly when neutrons are added: from 
2087 kev in ^^Cr to 1057 kev in "Cr, and down to 503 
kev in ^^Cr. The calculation reproduces this systemat- 
ics very well. This rapid decrease of the 9/2^ excita- 
tion energy can be understood as due to the attractive 
monopole interaction between the i'gg/2 and or- 



bitals that pushes down the vg^ /2 orbital. The high-spin 
states with positive-parity in 55,57(^j. g^j.g ^^^^ reproduced 
well. In both nuclei, the rotational bands built on the 
9/2^ state provide a strong experimental evidence for 
the shape-driving potential of this intruder level [l^ . 
In fact, we can see large E2 transitions in this band in 
later discussions. 

At this point, it is worthwhile discussing systematics 
of the 9/2f state. In Fig. [H the 9/2 states of the odd- 
mass Cr are compared with those of Fe, Ni, Zn, and Ge 
isotopes |31|]. As already seen from Fig. [5l the calcu- 
lated energy levels of 9/2^ are in a good agreement with 
experiment for ^^~^^Cr. Now in Fig. [SI all the curves 
show that the 9/2^ state energy steadily decreases with 
increasing neutron number. It should be pointed out that 
the more rapid drop of the 9/2^ energy from ^^Cr to ^^Cr 
correlates with the sudden increase of the calculated Oj 
energy from ^*Cr to ^°Cr in Fig. [3] Further calculations 
beyond = 37 predict a low-excited 9/2+ for ^i-escr. 
Interestingly enough, the curve for Cr isotopes is found 
to be very close to that of Ge isotopes, namely, the odd- 
mass Cr and Ge isotones have very similar 9/2^ energy. 
At A^ = 41, the 9/2+ state becomes the ground state for 
''^Ni and^^Qg^ .^^jj^g jg excited state at 0.158 MeV 
for ^^Zn. The calculation predicts that the 9/2^ level is 
an excited state for ^^Cr. 

In Fig. [3 we compare the experimental energy levels 
of 59 Cr andS°Cr with the present and the GXPFIA cal- 
culation. Our results in Fig. [71 as well as in the early 
discussed Fig. [31 show a much better agreement. The 
present calculation reproduces the positive-parity 9/2^ 
and 13/2^ states in ^^Cr, and predicts higher-spin states 
above them. For ^°Cr, our calculation nicely reproduces 
the experimental high-spin levels as well. It is thus clear 
that the vgg/2 excitations play a significant role in the 
current problem, and it is crucial to include this orbital 
for a correct description of these nuclei. 

The above success motivates us to move in on the even 
more neutron-rich region, where no experimental data are 
available. Let us next investigate the Cr isotopes beyond 
A^ = 38. Recent observation of very low first excited 2+ 
states in ^°Cr and ^^Cr suggests that these nuclei are de- 
formed 3. We thus expect that collective motion may 
be the dominant mode for the low-lying states in the iso- 
topes around A^ = 40. Figure [8] shows the calculated en- 
ergy levels of yrast states for the neutron-rich even-even 
isotopes 5^~^'*Cr. For those below A^ — 36, the calculated 
results reproduce well the observed trend of experimental 
data. In ^^Cr, the calculation correctly reproduces the 
enhanced excitation energy of the first excited 2^ and 4+ 
state as compared with the corresponding states in the 
neighboring ^''•^^Cr, which is, as discussed early, a strong 
signature of the subshell gap at A^ = 32. 

There is a striking behavior along the isotopic chain: 
the abrupt change in structure at A^ = 36 in ^°Cr. As can 
be clearly seen from Fig. [51 the 6^*", 8+, and 10+ state 
energies drop down suddenly at this neutron number. 
Another notable structure change, as shown in Fig. [HI is 
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FIG. 5; Experimental and calculated energy levels for the odd-mass ^^Cr. 
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FIG. 6: (Color online) Systematics of the first excited 9/2]^ 
states for Cr, Fe, Ni, Zn, and Ge isotopes. The calculated 
energy levels are compared with experimental data for odd- 
mass Cr nuclei. 



that the excitation energy of the state exhibits a peak 
at = 36. Remarkably, a similar structure evolution has 
been seen along the N = Z line in the proton-rich region, 
and a abrupt change in structure in the iV = Z = 36 nu- 
cleus ''^Kr has been discussed [1^ by the present authors. 
In Ref . , it was discussed that the yrast energy levels 
in the iV = Z < 36 nuclei ^"'Ge and ^®Se show larger sepa- 
rations, corresponding to smaller moments of inertia, but 
in ^^Kr, the level separation in the yrast sequence goes 
down drastically. We have suggested in Ref. [l^l that 
this is an example of phase transition along the N = Z 
line. 
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FIG. 7: Comparison of experimental energy levels with the 
present and the GXPFIA calculations for ^''•''°Cr. 



We now discuss the isotopes beyond N = 36. In the 
iV = 38 isotope ^^Cr, only the first excited 2f state has 
been observed [l^. Together with this state, a general 
trend starting from TV = 32 can be clearly seen, in which 



7 



8 - 



10^ 



> 
S 6h 



u 

0) 

=3 /I 

(U 4 - 

c 

o 



.6" 



24 



Cr 



8- 



cal 
exp 



28 30 32 34 36 38 40 42 44 46 



50 N 



FIG. 8: (color online) Experimental and calculated energy levels of the yrast states for the even-even Cr isotopes with TV 
28 - 50. 
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FIG. 9: (color online) Systematics of the first excited states 
for Cr, Fe, Ni, Zn, Ge, and Se isotopes. The calculated en- 
ergy levels are compared with the experimental data for Cr 
isotopes. 



the 2f excitation energies decrease gradually with in- 
creasing neutron number. The calculated excitation en- 
ergies of yrast states further drop down at the A'' = 38 
and 40 isotopes. The abrupt changes in structure across 
neutron-rich Cr isotopes appear as a result of low single- 
particle level density, which leads to a well-developed de- 
formation at = 40. When going toward more neutron- 
rich region from = 40 excitation energies of the yrast 
states begin to increase with increasing neutron number, 
which makes the TV = 40 isotope ^'^Cr the most deformed 
among the isotopic chain. It was proposed by Zuker et 



al. |22| that a minimal valence space to be able to de- 
velop quadrupole collectivity should contain at least a 
(j, J — 2, ....) sequence of orbits. In the upper (7g/2-shell 
nuclei, the 39/2 orbital and its quasi-SU(3) counterpart 
the orbital have to be taken into account in order 
to reproduce experimental data. Therefore, to describe 
the neutron upper f;g/2-shell nuclei, it would be better to 
take into account the orbital. At present, however, 
it is not possible to include the orbital to the shell 
model space because of huge dimensions of the configu- 
ration space. 



In Fig. [21 we show the calculated and experimental Oj 
states for even-even Cr isotopes, and for comparison, the 
experimentally known 0^ states for even-even Fe, Ni, Zn, 
Ge, and Se isotopes [31] are also plotted. For the only 
three Cr isotopes for which data exist, the calculation 
reproduces the data of ^''■^^Cr. The calculated Oj level 
at 5.6 MeV for ^^Cr (not shown in Fig. ^ lies too high 
when compared with the experimental one. The reason 
for this discrepancy may be attributed to the fact that 
neutrons in the /7/2 orbital are frozen in the calculation. 
From the systematics in Fig. [HI we can see two character- 
istics for these first excited Oj levels. One is the presence 
of a peak at TV = 36 in the experimental data for Ni and 
Zn isotopes as well as in the calculated results for Cr iso- 
topes. For Ge and Se isotopes, TV = 36 might also show a 
Oj peak; however, it is not conclusive at present because 
there are no data for the Oj energy at TV = 34. The other 
characteristic in Fig. [9]is the occurrence of a minimum in 
the Oj energy level around TV = 40. This minimum oc- 
curs because of the monopole interaction between i:f'j/2 
and 1^/5/2. 
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FIG. 10: (color online) Calculated and experimental values of 
(upper panel) B{E2,0i — » 2^), and (lower panel) i?(2f) for 
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FIG. 11: (color online) Calculated and experimental spectro- 
scopic qudrupole moment for even-even Cr isotopes. 



B. B{E2) values 

Calculation of electromagnetic transition probabilities 
serves as a strict test for wave functions of theoretical 
models. In the upper panel of Fig. IIOI experimental 
B{E2; 2+ 0+) dada for even-even Cr, Ni, Zn, and Ge 
isotopes are collected [III [3 IMI, and our shell-model 
calculations for Cr isotopes are compared with the known 
data. In lower panel of Fig. [TOl excitation energies of 
the first 2+ state E{2^) are plotted. With the usual 
effective charges = 1.5e and = 0.5e, the present 



calculation including the nf^ /2 orbital in the model space 
is found in good agreement with data for ^^■^''Cr. A sharp 
contrasting behavior between Ni and other isotopes can 
be clearly seen around TV = 40. The unusual B{E2) 
and E{2'^) behavior near ^^Ni shows up at the subshell 
closure TV = 40. The B{E2) values around = 40 are 
determined by filling the unique parity h'gg/2 orbital and 
the proton core polarization. The present calculation for 
Cr isotopes suggests that the predicted trend in B{E2) 
around TV = 40 is quite similar to the experimentally 
known systematics in Zn and Ge isotopes [ll|, [l2| . 

For shell model calculations in the model space con- 
sisting of the lp3/2, 0/5/2, lpi/2, 0^9/2 orbitals for pro- 
tons and neutrons with an inert ^^Ni core, large effective 
charges — 1.9e and — 0.9e were used to obtain 
the experimental B{E2) values for Zn and Ge isotopes 
pH [Ta [Tsj . This implies that excitation from the nf^ /2 
orbital would be important for obtaining large B{E2) 
values around TV = 40 Q. It should be noted that the 
calculated E2 transitions for 60,62,64qj. yjgjjj nearly the 
same values as the recent shell model calculations in the 
fpgd space including the ^5/2 orbital 14], where their 
B{E2) values are, respectively, 20.64, 20.7, 20.9 W.u. for 
60,62,64(^j.^ The calculated 2+ energies in Ref. [11] diverge 
from the experimental trend as they remain constant up 
to TV = 38 and increase at ^'*Cr. 

The deformation estimated from the calculated 
B{E2; 2+ ^ 0+) around TV = 40 is ^2 0.3, consis- 
tent with those obtained from the empirical formula [l3 ] 
and the Skyrme HFB calculations [lllli]. In the lower 
panel of Fig. [TOl we can see the TV = 32 gap in the 
Cr isotopes. With increasing neutron number beyond 
TV = 32, the first excited 2^ states decrease steadily to- 
ward TV = 40. The 2^ energy systematics of Cr behaves 
quite differently from Zn and Ge isotopes. The sudden 
drop in B{E2) values and large 2^ excitation energies for 
70-74qj. jj^g^y (j^^g -j-j^g neglect of the ^5/2 orbital [^ 
in the present model space. 

Figure [TT] shows a comparison for the calculated 
and experimental spectroscopic quadrupole moments Qs- 
The calculated Qs values for 52,54(^j, agree well with the 
data. The drastic variation in the calculated Qs values 
for ^""^''Cr may also be an artefact due to the neglect 
of the d^f2 orbital, which causes variations in the B{E2) 
values and 2^ energy systematics in these nuclei, as seen 
in Fig. \W\ 

In a previous paper [^ , we predicted a new band built 
on the first excited Oj state in two neutron-rich nuclei 
^^Ni and ^°Zn. The structure of states in this band is 
dominated by 2p-2h excitations from the fp shell to the 
intruder gg/2 orbital. We expect that an analogous band 
exists also in ^''Cr with TV — 40. The structure of this an- 
ticipated band is quite different from that of the ground 
state. This happens because the opposite signs of par- 
ity between the 179/2 orbital and the fp shell do not fa- 
vor Ip-lh excitations. Figure [T^] shows the theoretical 
level schemes. The predicted B{E2) values for ^^Cr and 
^*Cr are summarized in Table I. One can see strong en- 
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FIG. 13: Ground-state energies in the shell-model calculation 
for even-even Cr isotopes. The dotted curve is an extrapola- 
tion of the lower mass part of the results to show clearly the 
discontinuity in the slope of the curve at A'^ = 40. 



C. Ground-state energies and occupation numbers 



hancement in B(E2) values in the side band built on the 

0+) and 



first excited 0+ state in ^^-^Cr. The B{E2; 



B{E2;Al 



values for 



62,64 



2+ 
^1 



Cr are consistent with 



the calculated values by Sorlin et al. Q. This situa- 
tion is very similar to that in the neutron-rich ^^Ni and 
^'^Zn. The relevant ingredient here is excitations from 
the fp shell to the intruder orbital, and we expect 
that this is a common feature for the isotones of iV = 40. 
In Table II, B{E2) values for the positive-parity band 
built on the 9/2^ state in the odd-mass nuclei, ^^~^^Cr, 
are shown. We can see that for ^^Cr E2 transition rates 
have large values, and therefore show a strong collec- 
tivity. The deformation estimated from the calculated 
results is P2 — 0.18, consistent with the Total Routhian 
Surface (TRS) calculations 



TABLE I: Calculated B{E2) values for the positive-parity 
yrast states and the excited states in ^^Cr and ^^Cr. These 
values are compared with those of Sorlin et al. [sj. 
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It is striking that the effective SPE levels of the ^'P3/2 
and 1^/5/2 orbitals decrease drastically with increasing 
neutron number beyond TV = 42. This affects the prop- 
erties of the ground-state energy in the shell-model cal- 
culations. Indeed, as can be clearly seen in Fig. [131 the 
ground-state energies bend at iV = 40 so that it exhibits 
a discontinuity in the slope of the curve in Fig. 1131 This 
would correspond to a kind of phase transition at = 40. 

We can look at the above findings with other physi- 
cal quantities such as neutron occupation numbers of the 
fpg-sheW orbits for the low- lying 0^, 2+, and Oj levels. 
The upper panel in Fig. [14] shows occupation numbers 
for the ground states . For nuclei with smaller TV, 
the main component is of the VP3/2 orbital. At A'^ = 36 
and 38, the occupation number for the i'gg/2 orbital in- 
creases and nearly equals to that of the orbital. At 
A'^ = 40, the three occupation numbers of ^'^3/2, vgg/2i 
and J^/5/2 take almost same values. Going to higher neu- 
tron numbers from TV — 40, occupation number of the 
vgg /2 orbital increases drastically and becomes the dom- 
inant component. The occupation numbers for the first 
excited 2+ states in the middle panel show a similar dis- 
tribution as in the ground states. The bottom panel gives 
occupation numbers for the first excited 0^ states. At 



TABLE II: Calculated B{E2) values for the positive-parity 
band built on the 9/2+ state in the odd-mass nuclei 5^-5" Cr. 
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FIG. 14: (color online) Neutron occupation numbers of the 
/pg-shell orbits for three low-lying states in even-even Cr iso- 
topes. (Upper panel) Q'^ states; (middle panel) 2^ states; 
(lower panel) states. 



= 36, we can see a peak of the j^59/2 orbital in oc- 
cupation number, and the and J^59/2 orbitals take 
almost the same occupation numbers. This means that 
the Oj state in ™Cr consists of excitations from the i'Pz/2 
and vg%/2 orbitals, and lies higher than that of N — 34. 
For A'' = 38, the three orbitals provide the collectivity. 
For the Cr isotopes beyond N = 40, the orbital is 
the dominant component. 



IV. CONCLUSIONS 

We have carried out a comprehensive study for the 
structure of neutron-rich Cr isotopes using the spheri- 
cal shell model. Overall, the calculated results are in 



good agreement with the experimentally known energy 
levels for even-even 52-62(-;j. odd-mass 53-59Qj. jgQ_ 
topes. Furthermore, a lowering of excitation energies for 
neutron-rich Cr isotopes at and beyond TV = 40 has been 
predicted. We have explained several characteristic be- 
haviors found in the calculation: (1) a rapid decrease 
of the 9/2^ energies for the odd-mass Cr isotopes with 
= 29 — 35, (2) drastic structure changes at = 36 
with enhancement of the Oj energy and sudden reduction 
of the yrast band energy, and (3) steady decrease of the 
2^ energy in the smaller A^ region and the occurrence of 
pronounced collectivity around A'^ = 40. Several inter- 
esting consequences associated with the subshell closure 
A^ = 40 have been discussed: (1) enhanced B{E2) values, 
(2) lowering of the Oj energies, and (3) a new band with 
strong collectivity built on the Oj state in ^''Cr. 

We have pointed out that the key ingredients to pre- 
dict and explain these striking features are the monopole 
interactions. These attractive and repulsive interactions 
affect cooperatively the effective SPE. As we have seen in 
Fig. [21 the effective single-particle levels of 1^3/2, O/5/2, 
and 0(79/2 approach to each other near A^ = 40, due to the 
monopole effects. It seems that the cooperative effects of 
this behavior and the many-body correlations enhance 
the collectivity near A' = 40. On the other hand, the 
calculation without 1^5/2 (see Refs. 14| and [2^) does 
not lead to sufficient collectivity, which may suggest a 
need of including the 1^5/2 orbit in shell model calcula- 
tions. It is an interesting open question why the problem 
has not settled down and how one can unify the different 
views. The attractive 7^/7/2-1^15/2 lowers the f^/2 orbital 



with increasing A^, and also provides the mechanism of 
lowering the Oj energy around A^ = 40. The attrac- 
tive ^'^3/2-^^59/2 monopole interaction explains the rapid 
changes of structure at A^ = 36. Ge isotopes are consid- 
ered as the mirror nuclei of Cr isotopes with respect to 
the Z — 28 shell closure. We have pointed out that the 
systematics of energy levels and B{E2) values in Ge iso- 
topes is quite similar to that in neutron-rich Cr isotopes. 
We thus expect that our model can explain the charac- 
teristic behaviors in Ge isotopes, and are applicable to 
the N = Z nuclei with nuclear mass A = 68 — 100. 
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